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Tunneling and theoretical studies have suggested that Andreev bound states form at certain surfaces 
of unconventional superconductors. Through studies of the temperature and field dependence of the 
in-plane magnetic penetration depth A a & at low temperature, we have found strong evidence for 
the presence of these states in clean single crystal YBCO and BSCCO. Crystals cut to expose a 
[110] interface show a strong upturn in A a & at around 7 K, when the field is oriented so that the 
supercurrents flow around this surface. In YBCO this upturn is completely suppressed by a field of 
~ 0.1 T. 



PACS numbers: 74.25.Nf 

The low temperature behavior of the magnetic pene- 
tration depth (A) has played an important role in identi- 
fying the bulk order parameter in the high T c cuprates. 
Within the usual BCS quasiparticle picture, the temper- 
ature and field dependence of A are determined by the 
energy dependence of the quasiparticle density of states 
N(E), and hence give information about the magnitude 
of the energy gap, but not its phase. Recently however, it 
has been suggested that there is another contribution to 
A which arises from current carrying, zero energy, surface 
Andreev bound states (ABS) M, which form if there is 
ir phase shift between the different lobes of the order pa- 
rameter (2[jJ] . The observation of these states provides a 
key piece of evidence backing the d-wave scenario in the 
cuprates. In addition, this extra contribution has im- 
portant consequences for the interpretation of all surface 
impedance measurements in non-conventional supercon- 
ductors. 

The zero bias conductance peak (ZBCP) observed in 
tunneling measurements has been interpreted as result- 
ing from ABS Theoretically, ABS are expected to 
form most strongly on [110] surfaces, with the amplitude 
reducing monotonically to zero as the surface rotates to 
the [100], [010] or [001] directions. Experimentally how- 
ever ||, the ZBCP is seen for tunneling into all surfaces 
except [001]. This has been interpreted as resulting from 
surface roughness and/or nanofaceting [0. The contri- 
bution of the ABS to A can be understood simply by 
noting the increase in the penetration depth AA due to 
thermally excited quasiparticles is given by 



AA(T) 



N(E) df 
^ N(0) dE~ 



dE 



(1) 



where / is the Fermi function. For a ei-wave superconduc- 
tor N{E) ~ \E\, which gives AA(T)/A(0) = cvT, where a 
depends on the angular slope of the energy gap near the 
nodes Surface bound states add a singular contribu- 
tion 6(E) to N(E) which when substituted into Eq. dll) 



adds a divergent 1/T term to AA(T). The relative size of 
1 /T term depends on the orientation of the surface and 
the band structure of the material [0. The ABS contri- 
bution to A is highly non-linear in field, as the superflow 
Doppler shifts of the position of the bound state peaks. 
It is therefore expected that the 1 /T term is suppressed 
with a relatively weak magnetic field. 

In this letter, we report measurements of the temper- 
ature and field dependence of the in-plane penetration 
depth of high quality YBa 2 Cu306.g and Bi 2 Sr 2 CaCu208 
single crystals which provide strong evidence of the ex- 
istence of ABS in these materials. Evidence for an ABS 
contribution to A have previously been identified from 
studies of irradiated thin films H and grain boundary 
junctions 0. Here, we have been able to investigate the 
ABS more fully by measuring both the field and tem- 
perature dependence of A and by studying clean single 
crystals which give us more control over the direction of 
the Meissner currents and are free from intrinsic defects. 

We have identified the ABS contribution to A in YBCO 
by three different techniques. First, we have exploited the 
anticipated strong dependence of the ABS term on the 
crystal surface orientation by measuring several crystals 
which have been cut to expose different surfaces. Second, 
each crystal was measured in two orientations, with the 
probe field h parallel to either the afe-plane or the c-axis. 
Although the in-plane penetration depth is probed in 
both these orientations, the Meissner currents flow along 
quite different surfaces in the two cases. For h\\a or b the 
currents flow mostly along [001] surfaces where no ABS 
contribution is expected, whereas for h\\c the currents 
flow around the edges of the crystal and ABS contribu- 
tions are expected provided that a non [100]/[010] surface 
is present. This field orientation dependence also pro- 
vides an important check for impurity effects or isotropic 
paramagnetic contributions. Third, direct evidence for 
the presence of ABS states can be inferred from our mea- 
surements of the field dependence of \ a b- 
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Single crystals of YBCO were grown in yttria stabi- 
lized zirconia crucibles j^] , and were annealed for 3 weeks 
at 500° C in flowing oxygen to obtain optimal doping 
(T c ~ 94K, width ~ 0.2 K). Optical microscopy and x- 
ray diffraction showed that the crystals selected for this 
work were >90% untwinned. The BSCCO crystal was 
grown by the floating zone method and was annealed at 
350° C in 0.1% O2 to produce a slightly underdoped com- 
position with T c = 91 K ||]. Measurements of changes in 
A as a function of temperature and field were performed 
with resonant LC circuit driven by a tunnel diode oper- 
ating at 11.7 MHz [jn^, with frequency stability of a few 
parts in 10 10 Hz - 2. The dc field H is always collinear 
with the RF probe field h. 

The YBCO samples were measured in three geometries 
with the field applied parallel to either the a, b or c axes. 
For h parallel to either a or b the measured frequency 
shift is directly proportional to AA aj & [[l]J. In the third 
configuration (h\\c) the calibration is more difficult be- 
cause of the large demagnetizing factor. The calibration 
factor can be found either by normalizing the slope of 
the data to that for h\\a, b for 20 K< T <30 K where we 
expect the two configurations to be measuring the same 
quasiparticle response, or by using a recently developed 
approximate analytic solution for the screening currents 
in this geometry |lj] . We find that both methods agree to 
within experimental error (±5%). We denote the change 
in penetration depth relative to its value at our base tem- 
perature by AA^ b (T) and AX c ab (T) for h\\a,b and h\\c 
respectively. 

The temperature dependence of A a b for four samples 
of YBCO in the two field orientations is shown in Fig. 
|]. It can seen that for all samples AX^ b b (T) is highly 
linear at low temperature indicating high sample quality 
|^3|,0. However, the data in the second field orientation 
[AX c ab (T)] is markedly different below ~ 7 - 10K. There 
is also a clear correlation between the size of this upturn 
in A\^ b (T) and the shape of the sample. In the sequence 
of samples from A^D the proportion of non [100]/ [010] 
surface increases along with the size of upturn in A(T). 
These strong orientation dependencies of the upturn are 
clear evidence for it originating from ABS. 

Detailed calculations of the ABS contribution to A, in- 
cluding the effects of ABS broadening, have been per- 
formed by Barash et al. For a simple cylindrical 
Fermi surface the magnitude ((3) of the ABS 1/T term 
was calculated to be Uvfg(8)/[6kBX(0)], where g{8) = 
1 1 cos 3 (6) I — I sin 3 (8) \ | is the dependence of the upturn 
on the angle (8) the surface makes with the [110] direc- 
tion, and Vf is the Fermi velocity. The effect of a finite 
ABS lifetime is to broaden the (5-function contribution 
to N(E) making the divergence of A(T) weaker. The 
results of the numerical calculations in Ref. 0] can be 
approximated (for T 3> 0) by replacing the 1/T term by 
1/(T + T*). 




T[K] 

FIG. 1. Temperature dependence of the in-plane penetra- 
tion depths A\ a J(T) and AX c ab (T) for four crystals of YBCO. 
A constant has been added to AX^ b (T) so that is overlaps with 
AX^ b (T) at high temperature. The data for each crystal have 
been offset for clarity. The shape of the crystals is indicated 
in the figure. For sample D data taken down to 0.35 K in 
another apparatus is shown. The inset shows the difference 
between AA^ 6 (T) and A\% b (T) at 1.4 K versus the shape 
factor {g(8)y 

By subtracting the data for AX a ab (T) from AA^ b (T) 
we can isolate the ABS contribution. Fitting this to 
c/(T + T*) gives T* = 0.8±0.4K and c = 100 ± 30 AK 
for samples C and D, implying that the ABS states are 
only slightly broadened. The calculation of Barash et al. 
yields a somewhat larger value c = 10 3 AK for these crys- 
tals (assuming a cylindrical FS and nj=2x 10 5 ms _1 ), 
but we note that many factors, such as band-structure 
or diffuse scattering, can reduce c from this simple theo- 
retical estimate. 

We have attempted to test the correlation of the size of 
the upturn in A with the crystal shape in a more quantita- 
tive manner by calculating (g(8)) for our crystals, where 
() denotes an average around the surface of the crystal. 
8 was calculated from a digitized optical image of the 
crystal. The inset to Fig. shows AA^ fc (T) - AA^ 6 (T) 
at 1.4 K as a function of (g(8)). The linear fit shows 
that the correlation is good and the small intercept at 
(g(8))—0 indicates that there is not a significant effect of 
microscopic surface roughness below our optical resolu- 
tion. 

Next we turn to the field dependence of A. We have 
investigated this in two ways; either sweeping the tem- 
perature at fixed field or sweeping the field at fixed tem- 
perature. The results of the temperature sweeps at fixed 
field are shown in Fig. |^. A constant has been added to 
the AA° b (T) data in finite field so that they overlap with 
the zero field data at high temperature. By plotting the 
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FIG. 2. ZFC data for AX c ab {T,H) for sample D in fixed 
fields between and 10 mT, along with AA^ 6 (T) in zero field. 
The inset shows FC data for the same crystal. 



data in this way it can clearly be seen that an applied 
field suppresses the low temperature upturn in X(T). The 
good collapse of the data for T > 10 K implies that the 
field dependence of A [AX(H)], above 10 K is temperature 
independent. For each run the sample was cooled in zero 
field (ZFC) and then the field was applied at 1.4 K. The 
temperature was then cycled up to 30 K several times to 
check for drift and any hysteresis. For fields < 10 mT 
no hysteresis was seen and a zero field run following the 
finite field sweep overlap the virgin zero field curve ex- 
actly. We therefore conclude that for /iqH < 10 mT no 
vortices leak into the sample and we are in a reversible 
Meissner state. It can be seen that a field of ~ 10 mT is 
sufficient to suppress the upturn. For fields higher than 
~ 15 mT hysteresis was observed in the first T sweep 
after which the data became reversible. This hysteresis 
is presumably caused by flux initially leaking into the 
sample but becoming pinned away from the edges after 
the temperature cycle. Data were also taken where the 
sample was cooled through T c to low temperature with 
the field applied (denoted FC, inset Fig. ^). The effect 
of field in this case was much reduced (by approximately 
a factor of 4) . This is expected because of the more ho- 
mogenous flux profile and thus reduced surface field in 
this case. 

The effect of a dc magnetic field on the ABS contribu- 
tion can be understood as originating from a Doppler 
shift dE in the quasiparticle energy spectrum, where 
dE = evf ■ A (A is the vector potential). The effect 
of this shift on the ABS peaks has been observed di- 
rectly in tunneling measurements ^,|5|,[l5t. Substituting 
N(E) = S(E + ev f ■ A) into Eq. (@) gives 



A(0) 4T 



2k B T 



(2) 




showing that the 1/T divergence of A is suppressed when 
the applied field exceeds the temperature dependent field 
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FIG. 3. ZFC Field sweeps at fixed temperature for sample 
D. For two of the curves the 25 K data have been subtracted. 
The solid lines are fits to Eq. (El). 



scale H = k B T/(/j, evfX) = H T/T C , where H is of 
order the thermodynamic critical field [i§El c — <f>o/A£. 
A more exact formula for AX(H), which includes the 
necessary averaging over all quasiparticle trajectories has 
been given in Ref. Q. Numerically, the result does not 
differ significantly from Eq. (g) except in the effective 
value of H . There is also a contribution from the shift of 
the non-singular part of the quasiparticle spectrum but 
this is considerably weaker [AA/A(0) ~ H/Hq] [||. 

To make a detailed comparison to the data we first 
note that as the sample has a large aspect ratio, the ef- 
fective surface field will be significantly enhanced above 
the applied field. Also, we expect that the field will also 
vary somewhat across the thickness of the crystal, achiev- 
ing its maximum value near the top faces. We estimate 
an average 'demagnetizing' factor rj av — 1/(1 — N av ) by 
comparing the total frequency shifts when we extract the 
sample out of the coil in the H ||a, b and H\\c geometries. 
For sample D, r] av = 11.7. 

The results of field sweeps at fixed temperature are 
shown in Fig. [| At 1.34 K a very strong decrease in 
A with field is observed. Above 10K, AX(H) becomes 
positive and, to a good approximation, temperature in- 
dependent. For samples which do not show a strong 
upturn in AX^ b (T) (e.g., sample A) this behavior per- 
sists down to 1.3 K [|lCf| . We conclude therefore, that 
this behavior is a response of the bulk supercurrents. 
The response bears some resemblance to the prediction 
of Yip and Sauls for the non- linear Meissner effect [[|. 
However, the fact that it is largely temperature indepen- 
dent is in serious disagreement with the theory. This 
effect is nevertheless small in comparison to the total re- 
sponse at low temperature for samples C and D. The 
lower curve in Fig. || shows AA(H,T=1.35 K) with the 
temperature independent contribution subtracted along 
with a fit to Eq. (^. It can be seen that the agreement is 
excellent, with /io-ff = 1.9 mT. When allowance is made 
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FIG. 4. Field cooled measurements of AA„ b (T) for a 
BSCCO single crystal. The data have been offset for clar- 
ity. The inset shows the shape of the sample. 

for the field enhancement fioH e ff ~ 22.2 mT, implying 
Vf = 0.4 x 10 5 ms _1 . A fit to the Barash formula ||] gives 
a smaller value of fioH = 0.86 mT and a larger value of 
Vf = 0.86 x 10 5 ms _1 . At higher T we find that H in- 
creases but not as strongly as would be expected from Eq. 

. Qualitatively, this behavior was predicted by Barash 
et al. who showed that H tends to a constant as T — > T* . 
Thus at our lowest temperature H may be larger than it 
would be in the pure limit and the above analysis under- 
estimates the intrinsic value of Vf. With T* = 0.8 ±0.4 K 
the Barash formula gives Vf — 1.2 ± 0.2 x 10 ms . 

It is worth discussing briefly other processes which 
could effect the low temperature behavior of A. Small 
amounts of unitary scatters H], and non-local effects 
|fl6f can both cause \(T) ~ X^at low temperature, but 
do not cause an upturn. In some high T c cuprates a 1/T 
term in A(T) has been linked to the paramagnetic ions 
jl7j or amorphous surface layers jl|| . It is highly unlikely 
that such contributions would be present in clean single 
crystals of YBCO and in any case this effect would not 
depend on the orientation of the exposed surface or the 
field orientation in the manner observed here. Also the 
effect of field on such impurities would be much weaker 
than we have observed (e.g., ksT/g^B — 1 telsa, with 
3=2 and T=1A K). 

Measurements of BSCCO crystals show that similar 
effects occur in other cuprates. In Fig. || we show AA^ b 
for a sample of BSCCO in fields up to 97.8 mT. In zero 
field an upturn is clearly visible below ~ 1 K and again 
the application of held suppresses this upturn. Because of 
the very strong anisotropy of BSCCO it is not possible to 
measure X a b(T) in the H\\a,b configuration, although the 
absence of the upturn in this geometry precludes para- 
magnetic impurities. The very weak pinning in BSCCO 
also precludes ZFC measurements, as here we find ir- 
reversible behavior even in low fields. A larger field is 



required to suppress the upturn in this sample however 
this is largely explained by the smaller effect of a field 
in the FC case. This upturn has only been seen in this 
very clean sample of BSCCO. The smaller magnitude of 
the upturn [or the lower temperature of the minimum in 
A(T)] is explained by the surfaces of this crystal being 
close to [100]/[010]. 

In conclusion, by studying the effect of the orientation 
of the applied field to the crystal axes on the in-plane 
penetration depth, we have shown strong evidence for 
the existence of current carrying Andreev bound states 
on non [100] surfaces of YBCO. In clean crystal the states 
are only weakly broadened. The response is shown to be 
highly non linear in the Meissner phase, with a small 
field completely suppressing the anomaly. This behavior 
is consistent with a shift of the ABS peak in the quasipar- 
ticle density of states caused by the screening currents. 
Similar effects have been observed in BSCCO indicating 
that these effects are universal in cuprates and possibly 
all non-conventional superconductors. Besides being of 
intrinsic interest we also note that these effects must be 
taken into account in measurements of the surface prop- 
erties of all non-conventional superconductors which have 
surfaces where Andreev bound states can form. 

We acknowledge useful discussions with M. Aprili, D. 
Lee, J.R. Cooper and L.H. Greene. Work at Urbana was 
supported through State of Illinois - ICR funds. 
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